Enhancement of tunneling from a correlated 2D electron system by a many-electron 

Mossbauer-type recoil in a magnetic field 
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We consider the effect of electron correlations on tunneling from a 2D electron layer in a magnetic field 
parallel to the layer. A tunneling electron can exchange its momentum with other electrons, which leads 
to an exponential increase of the tunneling rate compared to the single-electron approximation. The effect 
depends on the interrelation between the dynamics of tunneling and momentum exchange. The results 
explain and provide a no parameter fit to the data on electrons on helium. We also discuss tunneling in 
semiconductor heterostructures. 



PACS numbers: 73.40. Gk, 73.21.-b, 73.50.Jt 

Low density two-dimensional electron systems (2DES) 
in semiconductor heterostructures and on liquid helium 
are among the most ideal many-electron systems. Such 
systems display strong effects of the electron-electron in- 
teraction, including those specifically related to electron 
correlations [jl]|^] . They show up dramatically in various 
unusual transport properties. One of the most broadly 
used techniques for investigating many-electron effects is 
tunneling a recent example being the observation Q 
of the giant increase of interlayer tunneling in double- 
layer heterostructures, apparently related to the onset of 
interlayer correlations. 

For electrons on helium, an exponentially strong devia- 
tion from the single-electron rate of tunneling transverse 
to a magnetic field has been known experimentally since 
1993 ||, but remained unexplained. Such a field couples 
the tunneling motion away from the 2DES to the in-plane 
degrees of freedom. The effect of the field and the role 
of electron correlations cannot be described by a simple 
phenomenological tunneling Hamiltonian. 

In this paper we provide a theory of tunneling from 
a correlated 2DES in a magnetic field B parallel to the 
electron layer. We show, using the model of a Wigner 
crystal (WC), that the tunneling is affected by the inter- 
electron momentum exchange and its dynamics, which 
is largely determined by short-range order. We discuss 
tunneling from 2DES on helium and in single quantum 
well heterostructures. The results explain and give a no 
parameter fit to the experimental data |)|, see Fig. 1. 
They suggest new types of experiments which involve 
tunneling through broad barriers and will be sensitive 
to short-range order in a 2DES. 

Electron correlations change the tunneling rate by ef- 
fectively decreasing the single-electron magnetic barrier. 
This barrier emerges because, when an electron tunnels 
from the layer (in the z-direction) , it acquires an in-plane 
Hall velocity vh = in the B x z direction and the 
corresponding in-plane kinetic energy mu>^z 2 /2, where 
lu c = eB/mc is the cyclotron frequency. Respectively, 



the energy for motion along the z-axis is decreased, or 
the tunneling barrier is increased by mwJz 2 /2. 

In a correlated 2DES, the tunneling electron exchanges 
its Hall momentum with other electrons, thus decreas- 
ing the energy loss This is somewhat similar to the 
Mossbauer effect where the momentum of a gamma quan- 
tum is given to the crystal as a whole Q] . In our case, the 
effect is very sensitive to the electron dynamics. If the 
rate of the interelectron momentum exchange u> p exceeds 
the reciprocal duration of underbarrier motion in imagi- 
nary time tJ 1 , then in- plane velocities of all electrons are 
nearly the same, and the Hall velocity is vjj cx 1/N — ► 
(N is the number of electrons). In this adiabatic limit 
the effect of the magnetic field on tunneling is fully com- 
pensated. For LOpTf ~ 1 a part of the tunneling energy 
goes to WC phonons, yet the B- induced suppression of 
tunneling is largely reduced. 
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FIG. 1. The rate of electron tunneling from helium surface 
W(B) as a function of the magnetic field B for the electron 
density n = 0.8 x 10 8 cm -2 and the calculated pulling field 
£± = 24.7 V/cm (solid curve). Lozenges show the experi- 
mental data |m. The error bars correspond to the uncertainty 
of the experimental parameters. The dotted curve is the cal- 
culation [^| for T = 0.04K without inter-electron momentum 
exchange. Inset: comparison of the present theory for B — 
with the experimentally measured density dependence of the 
tunneling rate. 
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In a strongly correlated system, where the electron 
wave functions overlap only weakly, one can "identify" 
the tunneling electron. Its out-of-plane motion for B = 
is described by the Hamiltonian 



Pz 

2m 



■U(z). 



(1) 



The potential U(z) has a well in which the electron oc- 
cupies the ground state, with energy E g . The well is 
separated by a tunneling barrier from extended states 
with a quasicontinuous spectrum, cf. Fig. 2 below. Wc 
assume that the tunneling length L is much less than the 
average inter-electron distance ~ n" 1 / 2 , where n is the 
electron density. Then small-amplitude in-plane electron 
vibrations about lattice sites are only weakly coupled to 
tunneling for B = We neglect this coupling. 

A magnetic field parallel to the electron layer mixes up 
the in-plane and out-of-plane motions. The Hamiltonian 
of the tunneling electron and phonons of the WC is H = 
H + H v + H B with 



and 



Hp 



2 2 

-muj„z 



■ OJrZN' 



-1/2 



Yfi x Pk ^ 



(2) 



(3) 



Here, pkj, Ukj, and lo^j are the momenta, displacements, 
and frequencies of the normal modes of the 2D Wigner 
crystal with the wave vector k, respectively (j = 1,2). 
We assumed that the equilibrium in-plane position of the 
tunneling electron is at the origin. Then its in-plane mo- 
mentum is p = N^ 1 / 2 Pkj- 

The Hamiltonian Hb couples the out-of-plane motion 
to lattice vibrations. The problem of many-electron tun- 
neling is thus mapped onto a familiar problem of a parti- 
cle coupled to a bath of harmonic oscillators ^Jul , with 
the coupling strength controlled by the magnetic field. 
However, there are two distinctions from the standard 
formulation. First, the coupling mixes together the par- 
ticle coordinate z and the momenta of the lattice. These 
two quantities have different symmetry with respect to 
time inversion. Because of broken time-reversal symme- 
try, the general problem of tunneling in a 3D potential in 
a magnetic field requires a special approach, which was 
developed earlier for an isolated particle [jy]. For the 
present model, the problem is simplified by the fact that 
in-plane motion is harmonic vibrations and the coupling 
is independent of |l(]] . 

The second distinction arises, because for 2DES the 
potential well U (z) is strongly nonparabolic near the min- 
imum (cf. Fig. 2). As a result, the standard instanton 
technique Jl2| does not apply |l3| ]. 

We will evaluate the tunneling rate in the WKB ap- 
proximation. In the presence of a magnetic field it is 



convenient to look for the WKB wave function under 
and behind the barrier in the momentum representation 
with respect to phonon variables, 



i>(z, {Pkj}) = exp[iS(z, {Pkj})], h=l, 



(4) 



and make a canonical transformation so that pkj and 
—Ukj be new canonical coordinates and momenta. 

To the lowest order in %, the action S in (||) can be 
obtained from the Hamiltonian equations for the trajec- 
tories of the system, 



S =p z z- y^UkjPkj, 
dH 



_ dH 

dpz' 

P^ = 



Pz 



dH 

dz 



dH 
du kj 



(5) 



In the (z, {pkj})-representation, the Hamiltonian equa- 
tions (^) have time-reversal symmetry. This allows us to 
solve them under the barrier in a standard way |l(| by 
keeping the coordinates z,pkj real and making the mo- 



menta p z , — Ukj, time t 



and action S(z, {Pkj}) 



ISe{z, {Pkj"}) purely imaginary. 

The Euclidean action Se(t) as a function of time is 
evaluated along a multidimensional trajectory (||) that 
goes under the barrier from the potential well to the 
boundary of the region which is classically allowed to 
both the tunneling electron and the WC vibrations. At 
this boundary one has to match the underbarrier solu- 
tion (with imaginary momenta) with the WKB solution 
behind the barrier (with real momenta), and therefore 



Pz(rf) = 0, u k j(r/) 



(6) 



where rj is the imaginary time at which the boundary is 
reached. 

We now discuss the initial conditions for the trajecto- 
ries (||). Typically, the characteristic intrawell localiza- 
tion length 1 /7 in the potential U(z) is small compared to 
the tunneling length L. For large jL ^> 1, the magnetic 
field may have strong cumulative effect on the tunneling 
rate, even where it only weakly perturbs the intrawell 
motion. Inside the well and close to it the electron in- 
plane and out-of-plane motions are then separated. We 
can set initial conditions at an arbitrary plane z — zq 
close to the well, yet deep enough under the barrier so 
that the wave function iJj(z, {pkj}) is semiclassical. For a 
harmonic WC, the dependence of ip on pkj is Gaussian. 
Then from (@) 



Se(0) = Y, p k j(0)p_ kj (0)/2ma; kj 

* — 'kj 



u k j(0) = -ip- k j(0)/mu> k j. 



(7) 



In the cases of interest, the dependence of ip on z is 
exponential near the well, ip oc exp(— jz). Therefore 



z(0)=z o , p z (0) = 17 = i^2m[U(z ) 



E n 



(8) 
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Under the barrier, the potential U(z) varies on the scale 
bigger than 1/7, and then 7 in (||) is independent of the 
exact position of the plane z = Zq. 

Solving the linear equations of motion ([5]) for the 
phonon variables Ukj , Pkj with the boundary conditions 
@> (0)i we can eliminate them, cf. Then Se takes 
the form of a retarded action for ID motion, 



U{z) = —Kz 1 - \e£±\z - muj 2 z 2 



(11) 



Se[z] = - 



2 Tf 



d,T\ 



^(dz/dr) 2 + U(z)-E g 



+-mu 2 z 2 (n) - (nuol/4N)J2^ 

kj 

x / rfr 2 z(Ti)z(T 2 )exp[-cj kj -(Ti - T 2 )] 
Jo 



(9) 



In (^|) we symmetrically continued the trajectory z(t) 
from Tf to 2t/, with z(rf + x) — z(rf — x) for < x < Tf, 
and set z = 0. The added section of the trajectory 
corresponds to underbarrier motion from the boundary of 
the classically accessible range back to the potential well. 
The tunneling rate W oc exp[— R], with R = 2minS , £. 

For small magnetic fields, the field-induced correction 
to the tunneling exponent (j^) is quadratic in u> c . It 
can be calculated along the zero-field trajectory dz/dr — 
[2/7(2) /m] 1 / 2 . This correction is always positive: mag- 
netic field decreases the tunneling rate. However, the 
correction is smaller than in the absence of the electron- 
electron interaction. 

Remarkably, although a part of the energy of the tun- 
neling electron goes to WC phonons, the tunneling rate 
increases with the increasing phonon frequencies. If the 
characteristic Wkj largely exceed the reciprocal tunneling 
time 1/r/, then z(t 2 ) ~ z(t%) in the second term in (^). 
As a result, the independent terms in (|)J) cancel each 
other, and tunneling is not affected by the magnetic field 
at all. This happens because, as the tunneling electron 
moves under the barrier, its in-plane momentum is adi- 
abatically transferred to the entire WC, similar to the 
Mossbauer effect. This can be contrasted with the case 
of an electron confined only inside the well but not un- 
der the barrier. Here the magnetic barrier is reduced by 
a factor of two compared to the free-electron case, but 
does not disappear ||. 

We now apply the results to electrons on helium and 
compare them with the experiment j^j. We will use the 
Einstein model of the WC in which all phonons have the 
same frequency u> p , which we set equal to the character- 
istic plasma frequency (27re 2 n 3//2 /m) 1 / 2 . The numerical 
results change only slightly when this frequency is var- 
ied within reasonable limits, e.g., is replaced by the root 
mean square frequency of the WC Q equal to |1J] 



1/2 



4A5e 2 n 3/2 /m 



1/2 



(10) 



For an electron which is pulled away from the helium 
surface by the field £±, the potential U(z) has the form 



for z > (outside the helium). On the helium surface 
(located at z = 0), U(z) has a high barrier ~ 1 eV which 
prevents the electron from penetrating into the helium. 

In (|ll|), the term oc A = e 2 (e - l)/4(e + 1) describes 
the image potential, e sa 1.057 is the dielectric constant. 
The field £± is determined by the helium cell geome- 
try and depends on the applied voltage and the electron 
density n, cf. |l6| . The term ex mui 2 = e 2 J2' |R-z | — 3 /2 
describes the Coulomb field created by other electrons 
at their lattice sites R; (the "correlation hole" |p|,^5[), 
for the tunneling length L < nT x l 2 . The conditions 
I/7 -C L -C vT 1 ! 2 are typically very well satisfied in 
the experiment, with 1/7 = 1/Ato 0.7 x 10 -6 cm, 
L = 7 2 /2m\e£±\ ~ 10~ 5 cm for typical £j_ ~ lOV/cm, 
and rT 1 ! 2 ~ 10~ 4 cm. 

The magnetic field dependence of the tunneling rate 
calculated from Eqs. (||) - is shown in Fig. 1. The 
actual calculation is largely simplified by the fact that, 
deep under the barrier, the image potential —A/z in ( |Tl| ) 
can be neglected. The equations of motion (||) become 
then linear, and the tunneling exponent R = 2SE(Tf) 
can be obtained in an explicit (although somewhat cum- 
bersome) form, which was used in Fig. 1. The correction 
to R from the image potential is ~ 1/7-k- When this 
and other corrections ~ 1/7-^ are taken into account, 
the theoretical curve in Fig. 1 slightly shifts down (by 
< 20% even for strong B), which is much less than the 
uncertainty in R due to the uncertainties in n and £±_ in 
the experiment ||. The theory is in excellent agreement 
with the experiment, with no adjustable parameters. 

The dependence of the potential U(z) on n gives rise 
to the density dependence of the escape rate W(B) even 
for B = 0. We calculated the exponent and the prefactor 
in W(0) by matching the WKB wave function under the 
barrier for 1/7 <C z <^ L with the intrawell solution. The 
latter was sought in the form t/j(z) = zexp[—A(z)]. The 
function dA/dz satisfies a Riccati equation which can be 
solved near the well (z <§; L) by considering the last two 
terms in ( |l~L| ) as a perturbation. When calculated to the 
first order in this perturbation, A allows to find not only 
the exponent, but also the leading term in the prefactor 
in the WKB wave function. The resulting tunneling rate 
is shown in the inset in Fig. 1. It fully agrees with the 
experiment [ p~7[ . 

For semiconductor heterostructures, tunneling in 
correlated systems has been investigated mostly for the 
magnetic field B perpendicular or nearly perpendicular 
to the electron layer, cf. The data on tunneling in 
a field parallel to the layer refer to high density 2DESs 
p8[ , where correlation effects are small. We expect that 
tunneling experiments on low-density 2DESs in parallel 
fields will reveal electron correlations not imposed by the 
magnetic field, give insight into electron dynamics, and 
possibly even reveal a transition from an electron fluid to 
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a pinned Wigner crystal with decreasing n. 

The effect of a parallel magnetic field is most pro- 
nounced in systems with shallow and broad barriers U (z). 
For example, in a GaAlAs structure with a square barrier 
of width L = 0.1 /xm and height 7 2 /2m = 0.02 eV, for 
the electron density n = 1.5 x 10 10 cm~ 2 and B = 1.2 T 
we have uj p tq « 0.6 and u c to w 1 (to = mL/"f is the 
tunneling duration for n = -B = 0). 
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FIG. 2. Relative rate of tunneling W = W(B)/W(0) vs 
magnetic field for a 2D WC in a semiconductor heterostruc- 
ture, with lutq = 0.5. Inset (a): the tunneling exponent R 
vs v — ^/2ujtq for lo c to = 1.0 (solid line) and 5 = (dashed 
line). Inset (b): the tunneling potential with (bold line) and 
without (thin line) barrier reduction due to static electron 
correlations. 

Electron correlations give rise to a coordinate- 
dependent lowering of the barrier, see Fig. 2. For nL 2 <C 
1, U(z) — 7 2 /2m — mu) 2 z 2 , < z < L [we count U 
off from the intrawell energy level E g ]. The picture of 
tunneling depends on the parameter v = \/2u>to. For 
v < 1 the electron comes out from the barrier at the 
point z = L where U (z) is discontinuous, cf. Fig. 2b. In 
this most important case, the boundary conditions (|^) 
for the tunneling trajectory should be changed to 



z ( T f) 



0. 



(12) 



but the tunneling exponent is still given by Eq. (|^). 

For B = the tunneling exponent R decreases with 
n, R = 7L[i/ _1 arcsin v + (1 — v 2 ) 1 ! 2 ] for v < 1, and 
R = tt^L/2v, for v > 1. Magnetic field causes R to 
increase and the tunneling rate to decrease. The effect is 
reduced by the inter-electron momentum exchange. The 
results for the Einstein model of the WC with uj-^j = u) p 
are shown in Fig. 2. The inset of Fig. 2 shows how R is 
decreased by the electron correlations even for B = 0. 

We have used the model of a WC to analyze the ef- 
fect of electron correlations on tunneling in a magnetic 
field parallel to the electron layer. We showed that the 
electron-electron interaction gives rise to an exponential 
increase of the tunneling rate compared to its single- 
electron value in a strong magnetic field. The effect is 
determined by the interrelation between the frequencies 
of in-plane electron vibrations and the reciprocal tunnel- 
ing time. For long tunneling time, the physics of large 



changes in the decay rate is closely tied to the physics of 
the recoilless fraction in the Mossbauer effect. Since the 
major contribution comes from the short-wavelength vi- 
brations, the results should apply not only to WCs, but 
also to all 2DESs with short-range order. Our results 
give a quantitative no-parameter fit to the experimental 
data H on tunneling of strongly correlated electrons on 
helium. 
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